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The Problem 


We introduce an unsolved combinatorial problem, which 
has been of recent interest to economists and political 
scientists. Then we present some new results, which were 
obtained from a computer simulation. 

The problem can be stated in the following scenario. A 
committee consisting of three judges A, B, and C wishes to 
choose, by a simple majority vote, one of three solutions 7 
or j or &. The preferences of each of the judges can be 
tabulated simply as: 


Judge A prefers 7 to j andj to k, 
Judge B prefers j to k and k to ¢, 
Judge C prefers k to 7 and 7 to j. 


If the committee votes on 7 and j, 7 will win; if the vote is 
on j and k, j will win; if the vote is on z and k, k will win. 
Thus the committee prefers ¢ to j, 7 to k, and k to 7. This 
committee can agree on no simple majority winner—what- 
ever issue is chosen, two of themembers prefer another to it. 
Such an impasse is called a cyclical majority. 

The unsolved combinatorial problem is: Given an odd 
number m of judges each with a transitive preference 
ordering over n (>3) issues, what is the probability that 
there will be no clear majority winner of the n issues when 
the committee votes on the issues in paired comparisons? 
Alternatively, how many of the (n!)” possible committee 
voting outcomes result in a cyclical majority? Three im- 
portant points should be kept in mind: 

(a) each judge has a strongly ordered transitive prefer- 
ence ordering over the issues—no judge is ever indifferent 
to any pair of distinct issues; 
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(b) each judge is equally likely to have one of the n! 
possible preference orderings over the 7 issues as his pref- 
erence ordering;’ 

(c) a cyclical majority is not said to have occurred when 
there is one issue that more than half the judges prefer to 
all others issues even though there may be cycles below 
first choice. For example, if in the preference pattern 
shown above there was an issue | that all three judges pre- 
ferred to 7, j, and k, then J would be a clear winner and 
the cyclical majority has not occurred in spite of the cycle 
within issues 7, 7, and k. 

This problem has a considerable history which we will 
not discuss here. The interested reader is referred to Black 
[2 (1948)], Arrow [1 (1951)], Riker [8 (1961)], Campbell 
and Tullock, [3 (1965)], and Klahr [6 (1966)]. We can 
summarize the relevant history by saying that many have 
tried to find a simple, closed-form expression for the prob- 
ability of a cyclical majority but none has succeeded.” 
Neither have we succeeded, but in our attempt we devised 
a method for simulating the voting process to estimate the 
probabilities which appears much more efficient than one 
previous simulation, Klahr [6], and more accurate than 
another, Campbell and Tullock [8]. 


The Results 


In Table I we show the fractions of a random sample of 
15,000 voting patterns that led to a cyclical majority for 
m judges and n issues, m = 3, 5, +--+ ,37,n = 3,4, +++ , 40. 
The fractions are given only to the nearest one-hundredth 
because, as Klahr [6] has shown, this is the theoretical 
accuracy limit for the probability when so small a sample is 
used. 

The same sample is used to generate all the entries in the 
table and these estimates are not independent of each other, 
as we explain below. Any one of the estimates of a cyclical 
majority from 15,000 patterns is accurate to more than .01 


1 This assumption is called ‘the assumption of an impartial cul- 
ture” by Garman and Kamien [5]. They also discuss the problem 
of partial cultures. 

2 It is, perhaps, an issue upon which reasonable men can agree to 
differ whether or not an expression such as > rer (simple expres- 
sion involving the vector r and the quantities m and 7) is asimple, 
closed-form expression, when to enumerate all ‘‘r € R” is currently 
an impossibly long computation for general m and n. Garman and 
Kamien [5], DeMeyer and Plott [4], and Niemi and Weisberg [7] 
provide formulae of varying complexity and use various combina- 
torial structures—primarily vectors, matrices, and multinomial 
coefficients—each unable to be evaluated for small m and n in any 
reasonable time. 
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PropaBiLity Estimates or a CyciicaL Masority. Tora SaMp.e, 15,000 Vorzs 


(Entries in table are not independent) 


TABLE 1. 
NUMBER 

OF 
ISSUES 3 5 7 9 ll 13 15 17 
3 205 207 07 207 98 208 208 08 
4 ell 014 014 216 216 216 o1l7 17 
5 01S 220 21 23 23 023 024 223 
6 20 025 27 028 28 029 229 230 
7 024 229 032 33 33 34 235 235 
8 227 033 36 38 38 238 39 39 
9 230 037 039 241 41 042 043 242 
10 232 239 243 045 045 045 046 246 
11 035 042 045 047 248 249 249 249 
12 236 245 248 049 250 251 251 051 
13 039 046 250 51 252 253 254 254 
14 o4t 048 052 53 054 255 255 256 
15 042 «50 254 255 056 257 258 258 
16 043 252 255 57 258 259 259 260 
17 045 253 57 59 260 60 261 61 
18 246 255 258 60 61 «62 062 63 
19 047 257 259 «62 062 263 063 264 
20 048 58 «6h 63 264 264% 265 265 
a1 049 59 062 064 265 265 266 266 
22 «50 260 263 265 266 266 267 067 
23 251 60 264 065 267 68 067 68 
24 51 261 065 266 268 269 269 269 
25 053 062 266 067 069 70 270 270 
26 054 063 «66 68 «69 70 71 71 
27 054 364 067 69 270 e771 Pres 72 
28 054 264 268 70 70 eT7l 12 73 
29 056 065 268 70 e771 072 073 273 
30 56 266 069 o71 072 73 73 73 
31 56 67 69 o71 072 73 0% 01% 
32 057 67 «70 272 73 07% 074% 274% 
33 58 267 70 272 73 74% 75 275 
34 58 ~68 eT 73 oT 075 275 16 
35 059 68 o72 074% 075 16 76 76 
36 059 269 72 74% 275 276 276 o7T 
37 260 «49 73 275 276 77 77 77 
38 260 70 073 275 «76 77 77 278 
39 061 70 14 076 tT «78 «78 78 
40 061 o71 014 16 77 278 78 18 


NUMBER OF JUDGES 


with probability. 95, but once this statement is made about 
one of the entries, the statistical accuracy of any other is 
not determinable without an involved conditional prob- 
ability analysis. In Tables II and III we provide two 
heuristic checks on the estimated cyclical majority prob- 
abilities. 

The exact probability calculated by Garman and Kamien 
[5] or by Sevcik [9] is shown to three decimal places in 
Table II.’ The paucity of numbers in Table II is some in- 


3 Garman and Kamien [5] and Sevcik [9] both counted the number 
of cycles, calculated the number (I!)”, and then calculated their 
ratio. Sevcik made the calculation for I = 3, 4, 5,6, 7 and J = 3, 
5, 7. Garman and Kamien made these same calculations (plus 
others) except for J = 6,7 andJ = 5,7. 
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dication of the difficulty of the problem. All the numbers 
for three issues are easily calculated with Garman and 
Kamien’s formula but, other than those shown, there are 
no precisely known numbers for issues greater than three. 
The sample 15,000 was generated in three separate 
groups of 5,000 each, and the probability estimates are 
available for each of the smaller samples. Some are shown 
in Table III. The smaller the variance among probability 
estimates over the replicated samples of 5,000 each, the 
more confidence there can be that the probabilities are 
generated by a stable process and that the numbers so 
obtained are good approximations to the true probabilities. 
The three sample estimates of any probability span a range 
less than or equal to .02 for 672 of the 684 numbers. For the 
twelve exceptions the range is .03. For example, the three 


Volume 13 / Number 4 / April, 1970 


TABLE II. Exact ProsBaBILity (To 3 PLACES) OF A 


CyciicaL Masoriry* 


Num- Number of Judges 

ber of 

Issues 3 5 7 9 ul eee 29 aie © 
3 .056 .069 .075 .078 .080 -085 -088 


4 411.1389 .150 6.156 = .160 
5.160 =.200)— 215 
6 .202 .251  .271 


7 .239 «6.296 = .318 


*From Garman and Kamien [5] and Seveik [9]. Numbers not 
shown are not known. 


estimates for 15 issues and 3 judges were .42, .43, and .40. 
In the data that generate Table I there are 105 entries 
identical in all three samples, 453 entries where two of the 
estimates are alike, and 114 entries where each of the three 
samples gives a different estimate. 

Our simulation provides several times as many cyclical 
majority probability estimates as any previous simulation. 
The cost of finding so many numbers is our inability to 
place statistical confidence limits on any but one of them. 
The reader can see that our estimates are never more than 
.01 away from the true probabilities when such are known; 
we think it reasonable to conclude that the others are, at 
worst, only slightly less accurate. Think of it this way: 
the estimate shown for 40 issues and 37 judges is known to 
be accurate within .01 with probability .95. The nature of 
the process is such that if there are systematic biases the 
greatest inaccuracy will be for 3 issues and 3 judges. But 
the benchmark comparisons indicate that the estimates for 
3 issues and 3 judges is also accurate to .01. 


The Method 


We now give our modus operandi for generating sample 
voting patterns and determining whether there is a majority 
winner or a cyclical majority. We reformulated the prob- 
lem for efficient computer simulation. 

For conceptual convenience we first describe the pro- 
cedure as it appears for a fixed number of issues and judges, 
and then we show how the evaluation is done simul- 
taneously for any lesser number ef judges or issues or both. 
Given a fixed number m of judges and fixed number n of 
issues, our computer program generates (15,000) ‘“out- 
comes,” defined as a random vote pattern for each judge 
and a count to determine the existence of a majority winner. 
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TABLE III. Ssiectep PROBABILITY CALCULATIONS OF 
A CycricaL Masority rrom THREE SAMPLES 
oF 5,000 
Num- Number of Judges 
ber of 
Issues 5 15 25 35 


10 .39, .39, .39 .47, .46, .46  .47, .47, 47 .49, .48, .47 
20 .59, .57, .57 .66, .64, .64 66, .65, .67 .67, .66, .65 


380 .67, .67, .64 .73, .74, .73 .74, .75, .76 .75, .75, .75 


40 .61, .60, .61 .77, .79, .79 .79, .80, .80 .80, .81, .79 


For each outcome, all judges but one are allowed to state 
their preferences at random. (The first judge always states 
his preferences in a fixed, canonical order—a mere labeling 
convenience. ) After all the preference orderings are known, 
each issue is tested as a possible winner. If no winner is 
found, a cycle is indicated. 

The preference orderings are developed in a preference 
vector. The testing for a winner is accomplished by develop- 
ing an overall preference matrix (first used by Garman and 
Kamien [5]) for each judge and then overall for the out- 
come. The entry in row r and column c of the matrix for a 
specific judge is 1 if issue ¢ is preferred to issue r and 0 
otherwise. The overall preference matrix is the sum of the 
individual preference matrices. The criterion used is that 
a majority winner exists if and only if there is a row in the 
overall preference matrix all of whose entries are less than 
m/2 (m = the number of judges). Two examples are shown 
in Figure 1, each with three judges (m = 3), A, B, and C, 
and four issues (n = 4), 7,7, k, and J. 

The preference order of the first judge, A, is fixed and the 
row and column labels for the overall preference matrix 
are chosen the same as in A’s preference matrix. As each 
judge’s preference vector becomes known, his individual 
preference matrix is generated and added to the overall 
preference matrix. After matrices have been generated for 
all judges and summed to form the overall matrix, each row 
of the overall matrix is scanned for an entry greater than 
m/2 (1.5 in the examples). If such an entry is found, the 
issue represented by that row cannot be a majority winner 
because there is some other issue that at least half the 
judges prefer. If some row has all its entries less than m/2, 
it is the majority winner’s row. Otherwise, each row has 
at least one entry greater than m/2; hence there is no clear 
majority winner and a cyclic majority has occurred. 

In the first example, issue & is the majority winner and 
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the overall preference matrix elements in the 7 and j rows 
which caused rejection of these rows are circled. In the 
second example, no majority winner exists and the ele- 
ments that forced each issue’s rejection are circled. 

Observation of our first efforts to program vote sampling 
showed that most of the computation time was spent in 
generating random preference vectors. For m judges and n 
issues, the vote generating process is proportional to mn’, 
while the evaluating process is proportional only to n’. 
We concluded that for a small fixed increment in processing 
time, the program could report, for each preference vector, 
not only the results for m judges and 7 issues, but also 
other results for any combination of judges up to m (odd) 
and issues up to 7. 

Recall that a preference matrix is developed for each 
judge and then a sum over all judges is computed. Before 
summing we have m preference matrices, one for each 
judge, over n issues. For any number m’ < m judges and 
n’ <n issues we can obtain an overall preference matrix 


(Ascending) preference 
order for each judge 1st Example 2nd Example 
over the issues 
Judge A t<j<k<il i<g<k<l 
B L<j<ic<k L<i<j<k 
Cc j<l<ick j<k<l<i 


_ Preference matrices for eachjudge: A“1”’ in row r and columne indicates that for this 
judge, the tssue represented in column c is preferred to that in row f. 


t j k t $ j k t 

A Oo 4 1 1 a as as | 
ji0 0 1 1 0 0 1 1 

ko 800 60061 0 0 0 1 

10 6 oO 0 0 0 0 0 

B i] 0 0 1 0 0 1 1 40 
j} 1 0 1 «=O 0 0 1 0 

ko 80 60 (OO 0 0 0 0 
jie 1 0 1 1 1 «0 

Cc io 0 1 0 0 0 0 0 
ji 1 0 1 21 1 01 1 

7 en i. 30.7 20: 7k 

wi 0 1 0 1 0 0 0 


Sum of preference matrices: (Circled elements eliminate that row issue as majority 
winner; row with * represents majority winner). 


| 0 1@ 1 0 @® 2 1 

J@ 0 38 2 1 0 @® 2 

*k| 0 0 0 1 1 0 0 ® 

ly] 2 1 2 0 ® 1 1 0 
issue k wins no winner 


Fig. 1. Two examples of program’s operation 
(83 judges, 4 issues) 
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using only the first m’ judges and the first n’ issues. The 
overall preference matrix constructed for m’ judges and n’ 
issues can itself be tested for a cyclical majority. Thus a 
given set of preferences for m judges over n issues can be 
used to generate a sample voting pattern for all possible 
judge-issue pairs (m’, n’) as long asm’ < mandn’ <n. 
The outcomes of the sample votes over various values of 
m’ and n’ are not independent of the outcome for m and n, 
but neither are they dependent in any known way. We 
believe it fair to say they are dependent in the same sense 
as outputs from a pseudo-random number generator are 
dependent. In general, a large body of data is reduced to 
one bit of information—a cyclic majority has or has not 
occurred. Reducing the number of judges or issues may 
change a cyclic majority into an outcome with a majority 
winner, or vice versa. 

The actual details of the machine language program for 
efficiently manipulating bits to make the proper counts are 
tedious and, for lack of general application, are not suit- 
able for presentation here. The program and details will 
be sent to interested readers on request. 
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